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Simultaneous  Nulling  in  the  Sum  ar  Oifference 
Patterns  of  a  Monopi  \ntenna 


1.  INTRODUCTION 

In  the  past  few  years,  considerable  research  and  development  has  been  expended 
on  adaptive  antennas.  Communications  and  sonar  systems  have  reaped  some  of  the 
benefits  of  adaptive  technology  while  radars  generally  lag  behind.  Some  of  the 
reasons  for  this  dichotomy  are  that  many  adaptive  techniques  are  not  suited  for 
microwave  frequencies;  that  phased  array  radars  are  large  with  many  antenna 
elements,  hence  more  adaptive  loops  are  required;  and  that  a  radar  has  tight  time 
constraints  for  detecting  and  tracking  targets  which  necessitates  the  use  of  fast 
response  adaptive  loops.  Consequently,  only  a  handful  of  radars  incorporating 
sidelobe  cancelling  techniques  exist  today. 

Monopulse  radars  add  additional  problems  in  the  implementation  of  adaptive 
processing.  The  problems  stem  from  the  coexistence  of  a  sum  and  difference 
pattern.  A  monopulse  phased  array  antenna  uses  two  antenna  patterns;  one  to 
detect  and  range  a  target  and  a  second  to  determine  a  target's  angular  location. 

Most  adaptive  antenna  research  has  not  directly  addressed  nulling  for  both  these 
patterns,  even  though  both  patterns  must  have  a  null  in  the  direction  of  interference 
to  enhance  the  radar's  performance.  Adjusting  the  far  field  sum  pattern  to  place  a 
null  in  the  direction  of  interference  will  not  place  a  null  in  the  difference  pattern  too. 
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Hence,  either  the  sum  and  difference  channeis  must  be  adapted  separately  or  a 
technique  of  simultaneously  nulling  in  both  channels  can  be  used.  Current 
approaches  to  this  problem  have  usually  assumed  the  first  approach,  and  the 
possibility  of  simultaneous  nulling  has  not  been  investigated. 

This  report  has  a  dual  purpose.  First,  it  describes  the  problem  of  nulling  in 
monopulse  antennas  and  points  out  that  both  the  sum  and  difference  channels  need 
to  be  considered  in  the  design  of  the  adaptive  processing  system.  And  second,  it 
theoretically  demonstrates  that  nulls  can  be  placed  simultaneously  in  the  sum  and 
difference  patterns  of  a  monopulse  phased  array  antenna  using  one  set  of  adaptive 
weights  shared  by  the  two  channels.  The  technique  used  to  generate  these  nulls 
offers  significant  hardware  and  software  savings  over  current  methods. 

The  next  section  of  this  report  introduces  the  concepts  used  in  array  analysis 
and  synthesis.  It  covers  linear  array  analysis,  sum  and  difference  patterns,  unit 
circle  representation,  and  low  sidelobe  synthesis.  This  section  provides  the  back¬ 
ground  and  notation  necessary  for  understanding  the  nulling  techniques  introduced 
later. 

The  following  two  sections  deal  with  nulling  in  sum  and  difference  channels. 
Section  3  performs  the  nulling  using  phase  and  amplitude  weights  while  Section  4 
shows  a  method  for  phase-only  nulling.  Both  techniques  are  demonstrated  for  the 
sum  channel,  then  the  difference  channel,  and  finally  simultaneous  nulling  in  both 
channels. 

Section  5  explains  the  computer  program  used  t_  demonstrate  the  above  analysis 
and  shows  several  results. 

The  final  section  outlines  the  pros  and  cons  of  using  the  techniq  :es  described. 

In  addition,  it  discusses  the  practical  implementation  of  this  method  and  future  work 
needed. 


2.  ARRAY  ANALYSIS  AND  SYNTHESIS 

This  section  describes  the  basic  concepts  required  to  understand  array  analysis 
and  synthesis.  The  theory  provides  the  basic  tools  needed  to  tackle  the  nulling  prob¬ 
lem  analyzed  in  the  following  sections.  A  great  deal  of  the  information  presented 

1  2 

here  was  drawn  from  array  antenna  books  by  Elliot  and  Steinberg.  ’  The  reader 
is  referred  to  these  books  for  additional  information. 


1,  Elliot,  Roberts.  (1981)  Antenna  Theory  and  Design,  Prentice-Hall,  Inc., 

Englewood  Cliffs,  N.J. 

2.  Steinberg,  B.  D.  (197  6)  Principles  of  Aperture  and  Array  System  Design, 

John  Wiley,  New  York. 
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Consider  a  linear  array  with  N  isotropic,  equally  spaced  (dQ)  radiating 
elements  (Figure  1).  Each  element  receives  signals  from  sources  in  its  environ¬ 
ment.  The  signals  pass  through  a  phase  shifter  and  amplitude  weight  to  a  summer. 
The  summer  adds  all  the  signals  together  to  generate  one  output  signal.  Adjusting 
the  phase  and  amplitude  of  the  received  signals  can  enhance  a  desired  signal  by 
rejecting  unwanted  signals. 


A  point  source  very  far  from  the  antenna  produces  a  uniform  plane  wave  across 
the  face  of  the  aperture.  If  the  source  is  at  broadside  (0  =  0°),  then  each  element 
sees  the  same  signal  phase  and  amplitude.  As  9  changes,  so  does  the  phase  seen 
by  each  element. 

The  difference  in  phase  between  two  elements  (m,n)  is 


A0 


m.  n 


kjf 


m, 


n 


(1) 


where 
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k  =  propagation  constant  =  2tt /  X 

A  =  wavelength 

l  =  the  additional  distance  the  uniform  plane  wave  must  travel 

m’  from  the  time  it  reaches  element  m  until  reaching  element  n. 

If  the  physical  center  of  the  array  is  also  the  phase  center,  then  the  phase  of  the 
signal  at  element  n  is  is 

4n  =  ^  dQ  [n  -  (N/2)  -  0.5]  sin  9  (2a) 

=  d  sin  9  =  d  u  (2b) 

n  n 


dn  =  kdQ  [n  -  (N/2)  -0.  5]  . 


(3) 


All  the  element  outputs  of  the  phased  array  are  summed  together  to  generate 
one  output.  At  broadside,  all  the  signals  are  equal,  so  the  output  is  equal  to  N 
times  the  signal  level  at  one  element.  If  the  signal  source  is  moved  in  an  arc 
around  the  array  (or  the  antenna  rotated)  the  output  of  the  summer  will  vary.  The 
phase  of  the  signal  at  each  element  changes  by  d  sin  9  as  9  varies.  Therefore, 
the  signals  do  not  add  in  phase.  In  fact,  as  9  increases,  a  point  is  reached  where 
all  the  signals  added  together  totally  cancel  and  result  in  a  null.  The  output  graphs 
an  interference  pattern  from  9  -  -90°  to  9  =  +90°.  This  far  field  antenna  pattern 
is  represented  by  the  summation  of  the  signal  amplitude  and  phase  at  each  element 


N 

S(u)  =  X 
n  =  1 


id  u 
J  n 


(4) 


The  mainbeam  lies  between  the  first  null  in  the  9  direction  and  the  first  null  in  the 
-0  direction.  Continuing  to  increase  9  produces  a  series  of  nulls  and  peaks.  The 
peaks  are  called  sidelobes  and  are  usually  much  smaller  than  the  mainbeam. 

The  maximum  output  from  the  summer  occurs  when  all  the  signals  entering  the 
summer  have  the  same  phase.  Normally,  this  condition  occurs  only  when  the  source 
is  at  broadside.  However,  the  signal  phase  at  each  element  may  be  changed  by  a 
phase  shifter.  Applying  a  phase  shift  of  -0  at  each  element  causes  all  the  signals 
to  have  the  same  phase  so  that  they  add  up  to  a  maximum  value  when  the  source  is 
at  some  angle  9  away  from  broadside.  In  effect,  the  linear  phase  shift  of  -C>n 
steers  the  mainbeam  in  the  direction  of  the  desired  source.  The  far  field  pattern 
of  the  steered  beam  is 


2 


(5) 


S(u)  = 


£  /dnU  -  V 
n=  1 


When  the  signals  at  each  element  are  added  in  phase,  a  sum  pattern  is  formed. 
A  mono  pulse  radar  uses  this  pattern  to  detect  and  determine  the  range  of  a  target. 
To  track  a  target,  a  monopulse  radar  uses  a  difference  pattern.  The  difference 
pattern  has  a  null  at  6  ■  0*.  It  is  formed  by  giving  the  signals  from  half  of  the 
array  elements  a  180*  phase  shift  and  adding  them  to  the  element  signals  from  the 
other  half  of  the  array.  The  resultant  pattern  has  a  different  sidelobe  and  null 
structure  than  the  sum  pattern.  Like  the  sum  beam,  the  difference  beam  can  be 
steered  to  a  given  direction  by  placing  a  linear  phase  shift,  -<f>  ,  across  the  array. 
The  far  field  difference  pattern  (N  even)  is  represented  by 
N/2  ....  N 


D(u) 


E  . 

n  =  1 


j(dnu  -< 


‘n>-  E 


j(dnu  -  0n) 


(6) 


N 

n-  -j  +  1 


The  sidelobe  level  of  the  sum  and  difference  patterns  are  controlled  by  weighting 
the  signal  amplitudes  at  each  element.  Equation  (5)  applies  to  a  uniformly  weighted 
aperture.  The  first  sidelobe  is  13.5  dB  below  the  peak  of  the  mainbeam.  A  differ¬ 
ence  pattern  with  a  uniform  amplitude  distribution  has  a  first  sidelobe  at  9.  5  dB 
below  the  peak  of  the  mainbeam.  Throughout  this  paper,  all  sidelobe  levels  will  be 
referenced  to  the  peak  power  of  the  far  field  array  pattern.  The  far  field  pattern 
of  a  weighted  aperture  is: 


S(u)  = 

N 

E 

a  /dnU 
n 

-♦n> 

* 

(7) 

n  =  1 

D(u)  * 

N 

E 

j(d  u 

b  e  n 
n 

-  V 

> 

(8) 

n  =  1 

where  an  =  sum  amplitude  weight  at  element  n,  and 

b  =  difference  amplitude  weight  at  element  n.  b  is  positive  when 
n  n  <  N/2  and  negative  when  n  >  N/2  +  1.  n 


The  sidelobe  reduction  is  not  without  cost.  Tapering  the  sum  amplitude  distribution 
lowers  the  antenna's  main  beam  gain. 
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In  most  cases  the  reduction  in  gain  is  small  and  the  lowering  of  the  sidelobes 
essential. 

A  useful  array  analytical  technique  is  the  unit  circle  representation  (Figure  2). 
The  circle  has  a  radius  of  one  and  is  centered  at  the  origin  of  the  real  and  imaginary 
axes.  This  circle  is  used  to  synthesize  a  far  field  antenna  pattern  in  the  same  way 
a  filter  response  is  designed  from  a  pole -zero  diagram. 


Equation  (5)  represents  the  far  field  sum  pattern.  Let 
jit' 

s  *  eJ^ 


(10) 


and 


U/  =  kd  sin  0  . 
o 

Substituting  these  values  into  the  equation  yields 
N 

S(u)  *  ^  a^  s(n  ^  (11) 

j.  aN-l  _N  -  2  !2a  +  li 


N  -  1 


Mil 


The  fundamental  theorem  of  algebra  says  a  polynomial  p(s)  with  N  -1  roots  can  be 
factored  into  the  form 


p(s)  ■  (s  -  s  j)  (s  -  s2>  (s  -  Sg)  . . .  (s  -  Sjj  _  j)  .  (13) 

The  polynomial  p(s)  represents  the  relative  far  field  pattern  of  an  array.  The 
N  -  1  roots,  Sj,  Sg.  ...  Sj^_  j,  are  the  locations  of  the  nulls  in  the  far  field  antenna 
pattern.  This  means  the  angular  location  of  the  nulls  is  given  by 

9n  =  sin-1  <*n/kdo)  .  (14) 

More  importantly,  the  coefficients  of  the  polynomial,  a^l a^  are  the  amplitude 
weights  at  each  element  in  the  array. 

The  angles  4/  determine  the  location  of  the  zeros  on  the  unit  circle.  This 
representation  gives  a  graphical  description  of  the  antenna  pattern  without  trans¬ 
forming  to  the  far  field.  Zeros  on  the  unit  circle  correspond  to  nulls  in  the  pattern. 
Approximately  half  way  between  two  zeros  is  a  lobe.  The  farther  the  zeros  are 
separated,  the  higher  the  lobe. 

Consider  a  uniformly  excited  array  with  four  elements.  The  far  field  pattern 
is 


.  .  3  2 

p(s)  =*  s  +  s  +  s 

+  1 

(15) 

(s+  1)  (s  +  j)  (s  -jl 

1 

(16) 

j*, 

■  *  =  -1 

4/1  =  n 

j'/'o 

•  =  -j 

=  3tt/2 

'  *  j 

*3  =  TTl  2 

As  long  as  the  coefficients  of  p(s)  are  real,  the  roots  occur  in  complex  conjugate 
pairs.  An  even  number  of  elements  always  has  a  root  at  \f/  ■  it.  Figures  3a  and  3b 
show  the  unit  circle  representation  for  this  array  and  its  corresponding  far  field 
pattern. 


AZIMUTH  AN9LE  IN  DEOREES 


Figure  3b.  Far  Field  Pattern  of  a  Four-Element  Uniform 
Array 


A  difference  pattern  may  be  analyzed  in  a  similar  fashion.  Consider  an  array 
with  the  far  field  pattern 

D(u)  =  pis)  =  s5  +  0.  41  s4  +  0.  59  s3  -  0.  59  s2  -  0.  41  s  -  1  (17) 

,  j37i/4.  .  -j37j/4.  .  -j 7T / 2 .  ,  j77/2.  .  i0. 

=  (s  -eJ  )  (s  -e  J  )  (s  -  e  J  '  )  (s  -eJ  '  )  (s  -  eJ  )  .  (18) 

The  zeros  of  the  pattern  occur  when  =  37i/4,  -3n/4,  tt/2,  -tt/2,  0.  Sidelobe  peaks 
are  approximately  at  7r/4,  -tt/4,  577/8,  -5)7/8.  A  diagram  of  the  unit  circle  repre¬ 
sentation  and  its  far  field  pattern  are  shown  in  Figures  4a  and  4b  respectively. 
Difference  patterns  always  have  a  zero  at  1 p  =  0". 


AZIMUTH  ANGLE  IN  DEGREES 


Figure  4b.  Far  Field  Pattern  of  a  Six-Element  Difference 
Array 


Although  the  past  two  examples  demonstrate  the  use  of  the  unit  circle  for  array 
analysis,  the  technique  is  generally  difficult  to  use.  Arrays  larger  than  four  ele¬ 
ments  present  problems  with  factoring  the  polynomial  p(s).  However,  the  unit 
circle  is  very  useful  in  array  synthesis. 

If  a  zero  indicates  a  null  location  in  the  pattern,  and  the  sidelobe  level  is  a 
function  of  the  distance  between  two  adjacent  zeros,  then  it  follows  that  by  maneuver¬ 
ing  the  zeros  on  the  unit  circle,  one  can  control  null  locations  and  sidelobe  levels. 
Moving  the  zeros  closer  together  in  one  area  separates  zeros  farther  apart  in  other 
areas.  To  make  all  the  sidelobes  lower,  all  the  zeros  are  moved  closer  to  the  nega¬ 
tive  real  axis  as  shown  in  Figure  5a.  Doing  this,  however,  increases  the  distance 
between  the  first  two  zeros  on  either  side  of  the  positive  real  axis.  These  zeros 
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determine  the  width  of  the  mainbeam.  Thus,  the  farther  the  zeroes  are  separated 
from  the  positive  real  axis,  the  lower  the  sidelobes  and  the  fatter  the  mainbeam. 


Im 


HIGH  SIDELOBES 
NARROW  MAINBEAM 


LOW  SIDELOBES 

WIDE  MAINBEAM 


Figure  5a.  Moving  the  Zeros  on  the  Unit  Circle  to  Get  Lower 
Sidelobes 


The  far  field  antenna  pattern  behaves  like  a  squeezed  balloon.  The  part  of  the 
balloon  that  is  squeezed  goes  down,  while  other  parts  of  the  balloon  expand. 
"Squeezing"  an  antenna  pattern  means  closely  spacing  the  zeros  on  a  unit  circle 
without  adding  any  new  zeros.  The  sidelobes  go  down  where  the  pattern  is  squeezed, 
but  other  parts  of  the  pattern  expand.  In  the  case  of  a  low  sidelobe  distribution  all 
the  sidelobes  are  "squeezed",  so  the  mainbeam  bulges. 

For  instance,  we  can  move  the  zeroes  of  the  four  element  uniform  array  closer 
to  the  negative  real  axis,  which  will  lower  the  pattern's  sidelobes  and  widen  the 
mainbeam.  Null  locations  at  \f/  =±  3it/4  instead  of  I (/  =  ±  n  12  give  a  far  field  pattern 
of 


S(u)  = 


(s  -  e^3ff/4>  <s  -  e'^4)  (s  +  l) 
s3  +  2.  4  s2  +  2.  4  s  +  1  . 


(19) 

(20) 


The  amplitude  weights  are  1,  2.4,  2.4,  and  1.  These  weights  give  the  far  field 
pattern  shown  in  Figure  5b.  This  diagram  confirms  our  prediction. 


AZIMUTH  ANGLE  IN  DEGREE8 


Figure  5b.  Far  Field  Pattern  of  a  Four-Element  Array  When 
the  Zeros  Are  Moved  on  the  Unit  Circle 


Placing  zeros  in  certain  predetermined  locations  can  control  the  sidelobes. 

Of  particular  interest  is  the  ability  to  specify  sidelobe  level  and/or  width  of  the 
mainbeam.  Some  important  low  sidelobe  sum  pattern  distributions  are  Chebychev, 

2  p 

Taylor,  Gaussian,  Cosine,  and  ( 1  -  r  )  .  The  Bayliss  distribution  is  the  most 
important  low  sidelobe  difference  pattern  distribution.  The  next  few  pages  describe 
how  to  generate  low  sidelobe  linear  array  distributions  using  the  unit  circle.  Only 
the  Chebychev,  Taylor,  and  Bayliss  distributions  are  considered. 

The  Chebychev  distribution  gives  sidelobes  of  equal  height  for  the  narrowest 
possible  mainbeam.  The  sidelobe  level  may  be  specified.  Chebychev  polynomials 
are  the  basis  for  the  distribution.  The  polynomials  are  characterized  by  equal 
ripples  of  a  specific  heignt.  Forcing  the  far  field  pattern  to  correspond  to  the 
polynomials  will  produce  the  desired  results.  The  Chebychev  polynomials  are  given 
by 


T^lu)  =  cos(kcos  1  v) 

-1  i  u  S  1 

(21) 

Tk<u)  =  cosh(kcos  *  v) 

1  u|  2  1  . 

(22) 

Making  the  variable  v  correspond  to  the  angle  u  will  relate  the  antenna  pattern 
S(u)  to  the  polynomial  T,  (v). 

1 

The  following  equations  perform  the  necessary  transformation. 

Tk(uQ>  =  b  (23) 
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< 


.  ..(Desired  sidelobe  level/ 20) 

b  =  10 

(24) 

u  =  cosh  |  1  [cosh  *  (b)]  > 

(25) 

°  |FT  f 

u^  =  cos  |(2r  -  1)  )j/[2(N  -  D]  j 

(26) 

i//  =  2  cos" 1  (u  /u  )  . 

“  r*  r* 

(27) 

When  the  number  of  elements  and  sidelobe  level  are  specified,  the  null  locations  on 
the  unit  circle  are  easily  determined.  Since  the  values  of  \j/  are  known,  the  poly¬ 
nomial  p(s)  in  factored  form  can  easily  be  determined.  Multiplying  all  the  terms 
together  gives  a  polynomial  of  degree  N-  1.  The  coefficients  of  the  s  terms  are  the 
amplitude  weights  for' the  array  elements. 

For  example,  consider  a  six-element  array  with  -20  dB  sidelobes. 

b  =  10(20/20>  =  10 

vq  =  cosh  (cosh  *  (10)]j  =  1.  1846 

vr  =  ±  0.  951.  ±  0.  588,  0 
l//  =  ±  73.  2°,  ±  120.  5°,  180°  . 

Next,  write  out  p(s)  in  factored  form. 
p(s)  =  ^s-e^73‘  2° 

Multiplying  these  factors  together  yields 

p(s)  =  s5  +  144s4  +  1.  85s3  +  1.  85s2  +  1. 44  s  +  1  .  (29) 

The  coefficients  of  p(s)  are  the  amplitude  weights  for  the  array  in  Figure  1.  This 
array  has  the  far  field  antenna  pattern  shown  in  Figure  6. 


)  (*-' 


J120.50 


)  H,81 


Figure  6.  Far  Field  Pattern  of  a  Six-Element  Linear 
Array  With  a  20  dB  Chebychev  Taper 


The  Chebychev  weighting  is  practical  only  for  a  small  linear  array.  As  N 

becomes  large,  the  amplitude  weightings  at  the  edge  of  the  aperture  become  large 

and  cause  mutual  coupling  problems  due  to  edge  effects.  Other  amplitude  distribu- 

3 

tions  are  better  suited  for  large  arrays. 

4 

One  such  distribution  was  developed  by  Taylor.  The  Taylor  distribution  was 
originally  developed  for  continuous  line  sources.  However,  the  technique  may  be 
extended  to  linear  and  planar  arrays.  For  a  linear  array,  the  line  source  is  sampled 
to  arrive  at  the  amplitude  weights  at  each  element.  A  technique  for  deriving  the 
element  weights  for  a  Taylor  distribution  is  explained  below. 

The  Taylor  distribution  is  similar  to  the  Chebychev  distribution  in  that  the 
maximum  sidelobe  level  can  be  specified.  The  difference  is  that  the  Taylor  distribu¬ 
tion  only  has  the  first  n  -  1  sidelobes  on  either  side  of  the  mainbeam  at  a  specified 
height.  All  remaining  sidelobes  decrease  exponentially  away  from  the  mainbeam. 

The  Taylor  distribution  may  be  thought  of  as  a  uniform  distribution  with  the  first 
n  -  1  zeros  on  the  unit  circle  moved.  The  low  sidelobes  are  obtained  by  moving  these 
zeros  closer  to  the  negative  real  axis.  Since  the  other  zeros  remain  untouched,  the 
outer  sidelobes  decrease  exponentially,  like  the  pattern  for  a  uniform  distribution. 
The  new  null  locations  are  given  by  the  formula* 


3.  Skolnik,  Merrill  L.  (1980)  Introduction  to  Radar  Systems,  McGraw-Hill  Book  Co., 

New  York. 

4.  Taylor,  T.  T.  (1955)  Design  of  line-source  antennas  for  narrow  beamwidth  and 

low  sidelobes,  IRE  Transactions-Antennas  and  Propagation,  AP^3: 16-28. 
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n 


(n  <  n) 


un  < 


j^A  +  (n  +  0.5)  J 


(n  £  n> 


(30) 


where  A  =  1/jr  cosh  *  (b).  These  null  locations  can  be  translated  to  null  locations 
on  the  unit  circle  by  the  formula 


^  =  2jr u  /N  .  (31) 

n  n 

From  tf/n,  the  factored  form  of  p(s)  is  formed.  Multiplying  the  terms  together 
gives  a  polynomial  whose  coefficients  are  the  Taylor  weights. 

A  20-element  linear  Taylor  distribution  with  n  =  5  and  -20  dB  sidelobes  has 
the  following  values 


u  =  ±  0.  55.  ±  1.  00.  ±  1.  89.  ±  3.  22,  ±  5,  ±  6,  ±  7  ±  8.  ±  9.  10 
n 

A  =  1.5  . 


Converting  ur  into  and  forming  the  polynomial  equation  yields 

p(s)  =  0.  67s19  +  0.  62s18  +  0.  59s17  +  0.  62s16  +  0.  72s15  +  0.  82s14 

+  0.  89s13  +  0. 93s12  +  0.  97s11  +  s10  +  s9  +  0.  97s8  +  0.  93s7  (32) 

+  0.  89s6  +  0.  82s5  +  0.  72s4  +  0.  62s3  +  0.  59s2  +  0.  62s  +  0.  67  . 

Figure  7  shows  the  resultant  far  field  pattern. 

A  difference  pattern  with  low  sidelobes  requires  a  new  set  of  weights.  The 
amplitude  distribution  for  sidelobes  in  a  sum  pattern  will  not  result  in  the  same  low 
sidelobes  in  a  difference  pattern  and  vice  versa.  For  instance,  when  the  Taylor  dis  - 
tribution  derived  above  is  used  in  the  difference  channel,  the  far  field  pattern  in  Figure  8 
is  obtained.  In  response  to  this  problem,  Bayliss  developeda  low  sidelobe  distribution 

5 

similar  to  the  Taylor  distribution.  The  Bayliss  distribution  is  characterized  by 
n-  1  equally  high  sidelobes  on  either  side  of  the  mainbeam  and  the  others  decrease 
exponentially  away  from  the  mainbeam, 

5.  Bayliss,  E.  T.  (1968)  Design  of  monopulse  antenna  difference  patterns  with  low 
sidelobes.  The  Bell  System  Technical  Journal.  May- June:  1968:623-650. 


Figure  T.  Far  Field  Pattern  of  a  20-Element  Array  With  a 
20  dB,  n  =  5  Taylor  Distribution 


1+0.  5>  2  6  _2 

A  +  n 


(n+  0.  5) 


,  2  2 
A  +  n 


A  +  n 


(1  S  n  £  4) 


(5  g  n  S  n  -  1) 


(n  £  n)  . 


The  values  for  6  are  given  in  Table  1.  The  variable  A  has  the  same  value  as 

it  did  for  the  Taylor  distribution.  is  found  by  using  Eq.  (31).  Finally,  the  poly 

nomial  coefficients  are  found  from  \b  . 

n 


Table  1.  Constants  for  Bayliss  Distribution J 


Bayliss 

Sidelobe  Level  in  dB 

Parameter 

15 

20 

25 

30 

35 

40 

A 

1.  0079 

1.  2247 

1.4355 

1.6413 

1. 8431 

2.0415 

61 

1.  5124 

1.6962 

1.  8826 

2. 0708 

2.  2602 

2.4504 

62 

2. 2561 

2. 3698 

2.4943 

2.6275 

2.7675 

2. 9123 

63 

3. 1693 

3. 2473 

3. 3351 

3.4314 

3.5352 

3.6452 

*4 

4. 1264 

4. 1854 

4. 2527 

4. 3276 

4.4093 

4.4973 

Consider  a  20-element  array  with  n  =  5  and  -20  dB  sidelobes.  The  values  for 


u  are 
n 


±  1.  81.  ±  2.  53,  ±  3.  47,  ±  4. 47,  ±5,  ±6,  ±7,  ±8,  ±9,  10  . 

In  turn,  the  polynomial  p(s)  is  found. 

p(s)  =  0.65s19+  0.71s18  +  0.82s17  +  93s16  +  s15  +  0.99s14 
+  0.90s13  +  0.72s12  +  0.46s11  +  0.  16s10  -  0.  16s9 

-  0.468  -  0.72s7  -  0.90s6  -  0.99s5  -  s4  -  0.93s3  -  0.82s2 

-  0.71s  -  0.65 

The  coefficients  of  p(s)  give  the  far  field  pattern  shown  in  Figure  9. 


AZIMUTH  ANGLE  IN  DEGREES 


Figure  9.  Far  Field  difference  Pattern  of  a  20-Element 
Array  With  a  20  dB,  n  =  5  Bayliss  Distribution 


The  remainder  of  this  report  describes  several  nulling  techniques.  The  unit 
circle  is  a  useful  tool  in  understanding  the  nulling  process.  An  antenna  has  only  a 
fixed  number  of  zeros  on  the  unit  circle  (called  degrees  of  freedom).  In  order  to 
place  a  null  in  a  desired  direction,  one  of  the  degrees  of  freedom  must  be  used. 

An  N  element  array  has  N  -  1  degrees  of  freedom.  Thus,  the  amplitude  weights  for 
the  20-element  Taylor  and  Bayliss  arrays  derived  earlier  have  only  19  degrees  of 
freedom. 


3.  AMPLITUDE  AND  PHASE  NULLING 

The  nulling  techniques  described  in  the  next  two  sections  are  based  on  the  beam 

g 

space  algorithm  described  by  Shore  and  Steyskal.  The  algorithm  is  modified  for 
difference  patterns  and  for  simultaneous  nulling  in  sum  and  difference  patterns.  The 
algorithm  is  one  that  minimally  perturbs  the  array  weights  (in  a  Least  Mean  Square 
sense)  when  forming  the  null.  In  this  way,  the  amplitude  and  phase  of  the  adaptive 
weights  produce  very  little  distortion  to  the  quiescent  far  field  pattern. 

The  algorithm  described  is  a  null  synthesis  technique  rather  than  a  closed  loop 
adaptive  algorithm.  The  location  of  the  jammers  must  be  known  for  the  nulls  to  be 
placed.  The  phase  and  amplitude  distribution  that  places  nulls  in  the  direction  of  the 
jammers  while  minimally  perturbing  the  adaptive  weights  is  then  computed.  The 


6.  Shore,  R.A.,  and  Steyskal,  H.  (1982)  Nulling  in  Linear  Array  Patterns  With 
Minimization  of  Weight  Perturbations,  RADC-  TR-82-32,  AD  A1 18695. 
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weights  used  to  generate  this  distribution  theoretically  yield  a  far  field  pattern  with 
the  desired  nulls.  In  a  practical  situation,  however,  antenna  errors  significantly 
limit  this  nulling  technique  unless  the  weights  are  adaptively  adjusted  using  some 
form  of  feedback.  It  is  possible  to  convert  this  open  loop  null  synthesis  algorithm 
into  a  closed  loop  adaptive  algorithm  to  compensate  for  the  antenna  errors.  How¬ 
ever,  that  practical  implementation  is  outside  the  scope  of  this  report. 

3.1  Sum  Channel 

Figure  1  serves  as  the  model  for  amplitude  and  phase  nulling  in  both  channels 
using  separate  weights  for  each  channel.  The  quiescent  weights  for  the  sum  channel 
are  given  by 

-j  d  u 

won  =  an  e  °  3  .  (35) 


Assume  that  a  weight  change  of  will  produce  the  required  nulls  in  the  far  field 
pattern.  The  new  cascaded  weights  are  represented  by 

w  =  w  (1+  Aw  )  .  (36) 

n  on  n 

The  weight  change  in  complex  form  is 


Aw  =  a  +  j  B 
n  n  J  H  n 


(37) 


When  there  are  no  weight  changes,  iw^  =  0  and  w^  =  wQn>  Substituting  Eqs.  (35) 
and  (37)  into  Eq.  (36)  leads  to 


-j  d  u 

w  =  a  e  ns(l+a+jfl). 
n  n  n  J  'n 


(38) 


This  new  amplitude  and  phase  distribution  puts  nulls  in  the  direction  of  the  jammers. 

The  far  field  sum  pattern  of  the  adapted  antenna  is  given  by 

N  .  .  .  . 

v- >  -jd  u  jd  u 

S(u)  =  2^ane  n  s  ( 1  +  a  +  jf}  )  e  n  (39) 

n=i  n  n  n 

N 

=  S  a_  e  “  °  +  2-,  aja_+j/3_)e  “  3  (40) 


jdn(u'us’  V  jdn(u'us> 

a  e  +  a  (a  +  i/3  )  e 

n= 1  n  n =  1  n  n  n 


N 

=  S  (u)  +  7.  a  (a  +  jfl  )  e 

n=  1 


j  d  (u  -  u  > 
J  n  s 


(41) 
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(42) 


N 

SQ(u)  =  Y  an  |cos  ldn<u*  us'5  +  3  ain  fdn(u‘  us^  I  ' 

n  =  I  ' 


and,  assuming  an  =  »N+  x  _  n  . 


N 

SQ(u)  =  Y  *n  cos  ldn<u  -  us>]  • 


(43) 


The  sine  terms  sum  to  zero  because  in  Eq.  (3a)  is  odd  symmetric  and  the  sine 
is  an  odd  function.  In  general,  M  nulls  are  required  in  the  pattern  to  cancel  M 
jammers  at  angular  locations  um  (1  g  m  =  M).  The  resultant  far  field  pattern  equals 
zero  when  u  =  u 


m 
N 

S(u  >  +  Y 

m  n  =  1 


a  (a  +  j  B  )  e 
n  n  n 


j  d  (u  -  u  ) 
J  n  m  s 


(441 


for  m  =  1,  2,  ....  M  . 

This  system  of  equations  may  be  looked  at  as  the  sum  of  the  quiescent  pattern  and 
M  cancellation  beams.  Each  cancellation  beam  matches  the  quiescent  pattern's 
amplitude  in  the  jammer  directions,  but  is  opposite  in  phase.  Therefore,  the  sum 
of  the  two  quantities  equals  zero  at  u^.  The  process  is  pictorially  represented 
in  Figure  10. 

The  above  system  of  equations  has  M  equations  and  N  unknowns.  To  determine 

the  values  of  Aw  ,  Eq.  (44)  is  solved  fora  and  /J  .  Shore  and  Steyskal  have  shown 
n  n  n  n 

E2  2 

(a  +  >3  )  minimizes,  in  the  least  mean  square 

n=  1  n 

sense,  the  adaptive  weight  perturbations.  This  is  a  very  desirable  characteristic, 
because  the  nulled  pattern  changes  very  little  from  its  quiescent  form. 

To  solve  the  system  of  equations,  we  put  it  into  matrix  form. 

AX  =  B  (45) 

where 
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SOLID  LINE  QUIESCENT  PATTERN 
DASHED  LINE  CANCELLATION  BEAN 


-OB  ■  M 


AZIMUTH  ANGLE  IN  DEGREES 

Figure  10.  Cancellation  Beam  and  Quiescent  Pattern 
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(A  A)  *  ( A  A )X  =  (A  A)  *  A  B 


X  =  Ar  (AAT)_1  B  , 


since  it  can  be  shown  that  (A^A)”1  A^  =  A^  (AA^)’1. 
The  product  of  AA^  is  a  matrix  whose  elements  are 


(AAf) 


jd  (u.  -u  )  jd  (u  -u  ) 

V  J  n  k  s  n  s  m 

a  e  a  e 

km  “  n  n 

n=  1 

V  2  Jdn<uk 


=  JC  4  |  cos(dn(uk  ~  um)l  +  J  sin  Idn<uk  "  um)^  | 

N 

=  z!  cos  |d  (u,  -  u  )]  for  k,  m  =  1,  2,  ....  M  . 
^ t  n  n  K  m 


The  sine  term  sums  to  zero  because  d  is  odd  symmetric,  the  sine  is  an  odd  func 

tion,  and  it  has  been  assumed  a  =  a., 

n  N  +  1  -  n 

Equation  (51)  may  be  written  as 

X  =  A^Y  (56a 


where 


Y  =  (AAV1  B  . 


The  vector  Y  is  found  by  inverting  (AA  )  and  post-multiplying  the  result  by  the 
vector  B.  Substituting  Y  into  Eq.  (51)  gives 
M 

«n+  J0n=  Z  yman  |cos!dn<us'VM  +  jsta[dn(ua-um)]|  *  ( 

m  =  1  ’ 


Next,  the  real  and  imaginary  parts  are  equated 
M 


a  =  A  y  a  cosld  (u  -  u  )1 
n  Jm  n  1  n  s  m  1 

m  =  1 


0  =  Z.  y  a  sinid  (u  -  u  )] 

n  •'m  n  1  n  s  m  1 


The  derived  weight  changes,  Awr  =  0<n  +  j$n»  are  used  to  adjust  the  complex  array 
weights  to  generate  the  nulls  in  the  pattern.  Using  Eq.  (36)  the  adapted  weights 
may  be  expressed  as 


w  =  w  y  expfjA  1 
n  on  'n  rlJrnJ 


where 


*n  =  l“n+1)2  +  ^l 


<t>n  =  tan*  [£„/<“„  +  Dl 


3.2  Difference  Channel 

The  nulling  technique  can  be  extended  to  include  difference  patterns.  Since  the 
sum  pattern  is  an  even  function  and  the  difference  pattern  an  odd  function,  the 
derivations  require  modification.  Weight  perturbations  that  create  a  desired  null 
in  the  sum  pattern  do  not  create  a  null  in  the  difference  pattern,  even  if  the  ampli¬ 
tude  distributions  are  uniform. 

One  way  of  producing  a  difference  pattern  is  shown  in  Figure  1  where  half  of 
the  element  signals  receive  a  180°  phase  shift.  The  quiescent  difference  channel 
weights  are 


-jd  u 
J  n  s 


w  =  b  e 
on  n 


As  before,  the  weight  change  Awn  produces  a  null  in  a  desired  direction. 


w  =  w  ( 1  +  Aw  ) 
n  on  n 


=  bne  ^\l  +  «n^n>. 


These  weights  have  a  far  field  pattern  represented  by 

A  idn(u"us>  A  *jdnUs  j  dn  u 

D<u)=2Jbne  n  S  +  2-bne  n8«*n  +  j0n>e  n 


V"  jd  (u-u  ) 

*Do<V* 

n=  l 


The  quiescent  far  field  pattern  is 


Dq(u)  =  ]C  |cos(dn(q- ug)]  +  j  sin[dn(u  -  ug)]  j 
n  =  1  " 


D(u)  is  zero  in  the  direction  of  the  jammers 
jd(u  -u  ) 

2]  (an+jfl)bne  n  m  s  =  -D  (u  ) 

„  ,  n  J  ^  n  n  o  m 


-  X,  i  b  cosfd  (u  -  u  )]  +  j  sin[d  (u  -  u  )]  j 
in  lnm  s'  J  1  n  m  sM 
n=  l  ' 


Because  b  is  an  odd  function 
n 


-  E 


b  cosld  (u  -  u  )]  =  0  . 
n  1  n  m  s  ‘ 


Equation  (68)  in  matrix  form  is  AX  =  B  where 


id2(ul_us) 


JdN(urus) 


jdl(u2-us) 


d2(u2  "  Us> 


JVVV 


.  p’idl(uM-us)  . 

bl  e  b2  e  ...  bN  e 


•>dN(uM’Us) 


B=  -[D^J.  Do(u2>.  ....  D^)]1 
X  =  [a  j  +  jp1#  a2  +  j02. 


The  method  of  least  squares  leads  to  a  solution 


X  =  A*  (AA^)  B 


=  A^  C 
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3.3  Simultaneous  Nulling  in  the  Sum  and  Difference  Channels 


Null  synthesis  for  sum  and  difference  patterns  in  the  manner  previously  des¬ 
cribed  requires  two  separate  sets  of  adaptive  weights.  Implementing  this  technique 
on  a  monopulse  radar  antenna  requires  N  complex  weights  in  the  sum  channel  and 
N  complex  weights  in  the  difference  channel  as  shown  in  Figure  1.  The  following 
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pages  describe  a  technique  that  allows  the  sum  and  difference  channels  to  share 
one  set  of  complex  weights  (Figure  11).  The  adaptive  weight  appears  before  the 
element  signal  is  split  into  the  sum  and  difference  channels.  Incorporating  such 
a  technique  into  an  adaptive  antenna  would  significantly  reduce  the  cost  and  com¬ 
plexity  of  the  system. 


Figure  11.  Monopulse  Linear  Array 


A  linear  array  with  N  elements  has  N  complex  weights  and  N-l  degrees  of 
freedom.  For  every  jammer,  the  sum  channel  uses  one  degree  of  freedom  and  the 
difference  channel  another.  The  drawback  to  this  technique  is  that  the  degrees  of 
freedom  available  for  nulling  in  the  sum  and  difference  channels  are  cut  in  half, 
so  that  (N-l)/2  degrees  of  freedom  are  available. 

Equations  (44)  and  (68)  still  hold  true,  but  now  they  must  be  solved  simultaneously. 

«  j  d  (U  ”  UQ)  V' 

/  .  (a  +  j  8  )  a  e  =  a  cos[(u  -  u  )]  (80) 

*  n  J  “  n  n  n  1  m  «  • 


a  cos[(u  -  u  )] 
n  1  m  s  1 


53  (0fn  +  JPn)bn  e  “ 


j  d  (u  -  u  ) 

J  r»  m  a 


-j  Z 


b  sin[(d  (u  -  u  )]  . 
n  1  n  m  s 


The  least  mean  square  solution  to  the  equation  is 


X  =  A  \  A  A  ^ )  B  . 


To  make  the  computations  easier  to  follow,  the  A  matrix  is  partitioned. 


A  = 

J 

K 

(AA^) 

J 

K 

JJ 

KJ 

[J  K  1 


The  elements  for  (AA  )  are 


=?4‘id",vv 


<AAf) 


i  =  a  b  e 

j  +  M.  k  “  n  n 


(AA 


' ) .  ,  ..  =  T]  a  b  e 

j,  k+  M  “  n  n 


j^j-v 


<AA  )  j+M,  k+M  —  -n 


t, 


Inverting  the  complex  matrix  (AA  )  and  multiplying  it  by  Y  results  in  2 
complex  matrix  R.  The  weight  perturbations  are 


an+^n=  AtR 


“jdjlu,  -  ug)  -jdl(uM-u) 

3 1  e  *1*  bl  e  ...  b,  e 


V 


(85) 

(86) 

(87) 

(88) 

(89) 

(90) 

(91) 

(92) 

2M  x  1 

(93) 


•jVrV 


■iV“ru.) 


-jdN<uM-ua) 


J  n  m 


M 


-j  d  (u  -  a  ) 

J  n  m  a 


-j  d  (u  -  u  ) 

*  n  m  a 


E,  ,  „  J  n  m  s  .  „  _  '  n  m  s  .  f 

(afe  +  1  a_  g_  e  +  b_t 

n  m 


m  =  1 


+  i  a  g  e 
J  n  6m 


n  m  +  M 


-jd^Um'V 

+  Jbngm+Me  > 


(95) 


(96) 


a  =  2^  la  f  cos(d  (u  -  u  )]  +  a  g  sin  [d  (u  -  u  >] 

n  .  |  n  m  1  n  m  s  n  “m  n  m  s 

m  =  l  ' 

+  b  f  ..  cosld  (u  -  u  )]  +  b  g  sinfd  (u  -  u  )  I 

n  m+M  1  n  m  sJ  n6m+M  1  n  m  s  | 

M 

B  =  2-i  I  -a  f  sinld  (u  -u  )]  +  a  cos[d  (u  -u  )] 

,  1  n  m  1  n  m  s  J  n  n  m  s’ 

m  =  l  ' 

-b  f  ..  sin[d  (u  -  u  )]  +  b  g  ..costd  (u  -u  )])  .  (97) 

n  m+M  1  n  m  s  1  n  m+M  1  n  m  s 


All  three  techniques  were  modeled  on  the  computer  and  the  results  are  shown 
in  Section  5. 


4.  PHASE-ONLY  NULLING 

Nulling  with  phase  shifters  is  much  more  desirable  than  nulling  with  complex 
weights.  The  hardware  for  phase-only  nulling  is  less  complicated.  Adjusting  the 
signal  amplitude  requires  an  attenuator  or  two  phase  shifters  in  a  bridge  circuit 
for  each  element.  The  extra  losses,  mismatch,  and  cost  of  the  amplitude  control 
discourages  its  use.  Also,  phased  arrays  usually  have  only  phase  shifters  and 
incorporation  of  amplitude  weighting  requires  modifications  to  existing  antennas. 
An  adaptive  technique  that  uses  phase  control  rather  than  phase  and  amplitude 
control  can  more  easily  be  implemented  on  existing  phased  arrays. 

4. 1  Sum  Channel 

Phase  only  nulling  in  the  sum  channel  of  a  low  sidelobe  antenna  is  discussed 

O 

in  the  literature.  The  technique  closely  follows  the  process  used  in  phase  and 
amplitude  nulling.  The  quiescent  sum  channel  weights  are  given  by  Eq.  (35).  To 
put  nulls  in  the  antenna  pattern,  a  phase  shift  of  0^  is  given  to  each  element. 

8.  Baird,  C.  A. ,  and  Kassweiler,  Ci.  C.  (  1 976)  Adaptive  nulling  using  digitally 
controlled  phase-shifters,  IEE  Trans,  on  Antennas  and  Propagation, 

24( No.  51:638-649. 
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a* 
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«r+ 
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(98) 


-id  u  10 
^  J  n  s  Jrn 
w  =  a  e  e 

n  n 


-  an  e  *****  (1  +j0n).  0n  «  1 


(99) 


The  approximation e  n  —  1  +  j<j>  is  necessary  to  solve  the  problem,  but  it  makes 
the  answer  approximate.  Fortunately,  this  approximation  works  very  well  for 

g 

low  sidelobe  antennas. 

The  far  field  pattern  of  these  weights  is  calculated  from 
N 


E-i d  u 
an  e  n  S 

n=  1 


(1+  j0n)  e 


i  d  u 
J  n 


N  .  .  /  v  N 

V'  J<Uu-u  J 


n=  1 


a  e 
n 


+  j  £  a  0  eJ  n 

J  “  n  rn 


idJu-u) 


n=  1 


(100) 


(101) 


At  the  jammer  angles,  the  pattern  goes  to  zero 


N  ...  .  N 

l  d  (u  -  u  ) 
a  \  V  n  m  s 

S(u  )  =  7,  a  e  + 

m  n  -  , 

n=  1  n=  1 


i d  (u  -  u  ) 
J  "  m  s 


V  *  J  n 
j  L  an  e 

n=  1 

1  £an«n  jcos|dn(Um-us)]+jsinldn<um-us>l  | 

n  =  1  ’ 

N 

=  -  7,  a  |  cos(d  (u  -  u  >]  -+  j  sin  Id  (u  -u  )]  l 

n~“j  n)  lnm  s'  J  1  rv  m  s  j 


(102) 


(103) 


Next,  equating  the  real  and  imaginary  parts,  assuming  an  real 
N  N 

Ea  0  sin  Id  (u  -u  )  =  7.  a  cos  Id  (u  -  u  )] 

r\  1  n  m  ft  n  n  m  S 


.  n  n  n  m  s  n 

n  =  1  n  =  1 


N 


N 


Ea  0  cos  d  u  -  u  )]  =  7.  a  sin  d  (u  -  u  )]  . 

.  n  n  1  n  m  s  1  •—  n  1  n  m  s 

n=  l  n=  1 


(104) 


(105) 


The  second  equation  equals  zero  because  the  sine  is  an  odd  function,  dn  has  odd 

symmetry,  and  a..  ,  =  a  has  been  assumed. 

N  +  1  n 

Equation  (104)  can  be  put  into  the  matrix  form  AX  =  B  where 


A  = 


aj  sinldjluj  -  us>]  ...  aN  sinid^Uj  -  ug)] 


a  j  sinldjlu^j  -  ug>  ...  aN  sinldN(uM  "  “s’1 


(106) 


37 


|yC  an  costdn(uM'Us)1J 

The  unknown  phase  shifts  are  obtained  from  the  method  of  least  squares. 


X  = 

T 

aac 

(109) 

c  = 

t  -  1 

(AAA  )  B 

(110) 

<AAT)jk  -  sin[(dn(u.-ug»  sin[dn(uk-us)] 

M 

0  =  2-i  c  a  sin[(d  (u  -  u  )]  . 

n  * — ,  m  n  1  n  m  s 

m  -  I 

The  variable  c  is  the  mth  element  in  the  vector  C.  Nulls  are  made  in  the  sum 
m 

pattern  in  the  directions  of  the  jammers  when  0 n  is  applied  to  the  phase  sh;fters. 

4.2  Difference  Channel 

The  derivation  of  the  phase  shifts  for  the  difference  channel  follows  that  of  the 
sum  channel. 

(113) 


-jd„u=  J' 


.  ^  n  s  n 

w  =  d  e  e 


“  bn  e 


-id  u 
J  n  s 


(111) 

(112) 


(1  *  J0  )  . 


The  far  field  difference  pattern  is 

*  j  d  (u  -  u  ) 

D(u)  =2^  b  <1  +  j0  >  e  n  s 
n=  1  n 


D(u)  equals  zero  when  u  =  u  . 

m 

N 

i  b  0  |  cosld  (u  -  u  )]  +  j  sinld  (u  -  u  )j } 

n  =  1  ' 

N 

=  -  /_  b  |  cosld  (u  -  u  ))  +  j  sinld  (u  -  u  )]  \ 
n  j  1  n  m  s‘  J  lnm  s  j 

Equating  the  real  and  imaginary  parts  of  the  equation, 

N  N 

Sb  0  cosld  (u  -u  ))  =  -  y.  b  sinld  (u  -  u  )) 

.  n  n  1  n  m  s'  x—‘.  n  1  n  m  s' 

n  =  1  n  =  1 


Eb  0  sinld  (u  -  u  )]  =  -  b  cosld  (u  -u  ))  ,  (117) 

,  n  n  L  n  m  s  1  ,  n  1  n  m  s' 


Equation  (116)  equals  zero. 

The  remaining  system  of  equations  can  be  solved  using  the  method  of  least 
squares. 


AX  =  B 


T  T  -1 
X  =  AX(AA  )  B 


-  atc 


(119a) 


(119b) 


bj  cos[dj(Ul  -ug)j  ...  bNcos[dN(Uj -ug)] 


bf  costdj^ -ug)]  ...  bNcos[dN(uM-us)]. 


-  y.  b  sinld  (u.  -u  )] 
n  1  n  1  s' 


Eb  sinfd  (u.  -  u  ) 
,  n  1  n'  M  s 
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B  = 


N 


T  ancos[dn(Ul-us)] 
n  =  1 


N 


£  a„cosl(VuM‘Us)l 

n  =  1 
N 

-  L  bnsin[dn(urus)l 


n  =  1 


N 


-  53  bnsinldn(uM-Us)] 

n  =  1 


The  least  squares  solution  is 


T 

X  =  AlC 


T  - 1 

C  =  (AA  )  B 
N 


(AA7)^  .  =  £  a^  sin(dn(u.  -  ug))  sin(dn(u .  -  ug)].  i.  j  £  M 

N 

(AAT)i  +  M,j  =  ^3  an  bn  sintdn(ui  'V1  COSldn(uj  'V’ 

n  =  i 

N 

(AAT)i  :  +  N  »  Z  anbncos[dn(u.-us)]sin(dn(u  -ug)] 
j  n  ~  1 

N 

<AAT)i  +  M.j  +  M  =  bn  cos[dn(u.-us)]  cos[dn(u.-us)] 

M  2M 

<p  -  ^  c  a  sin[d  (u  -  u  )J  53  c  b  cos[d  (u  -u  )] 

"  m  =  1  nmsm  =  MHmn  nms 


(125) 


(126) 

(127) 

(128) 

(129) 

(130) 

(131) 

(132) 


The  variables  c  are  elements  of  the  matrix  C.  These  values  of  <p  simultaneously 
m 

generate  a  null  in  the  sum  and  difference  patterns  in  the  direction  of  interference. 
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5.  SIMULATION  RESULTS 


I  modeled  the  nulling  techniques  described  in  Sections  2  through  4  on  the  com¬ 
puter.  Appendix  A  contains  a  block  diagram  and  description  of  the  program. 

This  section  shows  the  results  of  several  computer  runs  of  the  model.  The  nulling 
results  from  this  section  are  displayed  in  Table  2. 

The  antenna  array  was  modeled  after  Figure  11.  Both  the  sum  and  difference 
channels  share  one  set  of  adaptive  weights.  The  output  for  the  channel  is  obtained 
by  multiplying  the  received  signal  at  each  element  by  the  adapted  weight.  Next,  a 
3  dB  coupler  sends  half  the  signal  to  the  sum  channel  and  half  to  the  difference 
channel.  Each  sum  channel  signal  is  multiplied  by  an  amplitude  weight  then  summed 
to  obtain  one  output.  The  difference  signal  is  obtained  in  a  similar  manner,  except 
that  half  of  the  signals  are  phase  shifted  by  180°  before  being  summed. 

The  first  series  of  runs  used  a  10-element  uniformly  weighted  linear  array. 
Figures  12a  and  12b  show  the  quiescent  far  field  sum  and  difference  patterns  and 
weights  for  the  array.  Initially,  the  adaptive  amplitude  and  phase  weights  were  set 
to  1.  0.  When  a  jammer  was  placed  at  33°,  these  weights  changed  to  produce  a  null 
in  the  pattern  in  that  direction.  If  a  null  is  placed  only  jn  the  sum  channel,  the  re¬ 
sults  in  Figures  13a- 13d  are  obtained.  The  sum  channel  weights  in  Figure  13a  are 
obtained  by  multiplying  the  adapted  weights  by  the  sum  channel  amplitude  distribution. 


E 


Figure  12b.  Quiescent  Far  Field  Difference -Channel 
Pattern  of  a  10 -Element  Uniform  Array 


VARIABLE 

COMPLEX  WEIGHTS  AMPLITUDE  TAPERS 


ELEMENT 

AMPLITUDE 

PHASE 

SUM  CHAN 

DIF  CHANNEL 

1 

1.00000 

0.00000 

1.00000 

-1.00000 

2 

1.00000 

0.00000 

1.00000 

-1.00000 

3 

1.00000 

0.00000 

1.00000 

-1.00000 

4 

1.00000 

0.00000 

1.00000 

-1.00000 

5 

1.00000  ' 

0.00000 

1.00000 

-1.00000 

6 

1.00000 

0.00000 

1.00000 

1.00000 

7 

1.00000 

0.00000 

1.00000 

1.00000 

8 

1.00000 

0.00000 

1.00000 

1.00000 

9 

1.00000 

0.00000 

1.00000 

1.00000 

10 

1.00000 

0.00000 

1.00000 

1.00000 

Figure  12c.  Quiescent  Uniform  Sum  and  Difference  Channel  Weights 


The  difference  channel  weights  are  the  result  of  multiplying  the  adapted  amplitude 
weights  by  the  difference  channel  amplitude  distribution.  Adapted  phase  values 
are  the  same  for  both  channels.  Corresponding  to  the  new  adapted  weights  is  the 
far  field  sum  pattern  shown  in  Figure  13b.  Its  quiescent  pattern  and  associated 
cancellation  beam  are  shown  in  Figure  13c.  These  same  adapted  weights  do  not 
yield  a  null  in  the  difference  pattern  at  33"  (Figure  13d). 


Table  2.  Nulling  Results  for  Figures 


VARIABLE 

COMPLEX  WEIGHTS 

VARIABLE  WEIGHTS  TIMES 
AMPLITUDE  TAPERS 

ELEMENT 

AMPLITUDE 

PHASE 

SUM  CHAN 

DIF  CHANNEL 

1 

0.98951 

-0.10124 

0.98951 

-0.98951 

2 

0.90362 

0.03252 

0.90362 

-0.90362 

3 

1.04667 

0.08782 

1.04667 

-1.04667 

4 

1.08633 

-0.05069 

1.08633 

-1.08633 

5 

0.93674 

-0.08166 

0.93674 

-0.93674 

6 

0.93674 

0.08166 

0.93674 

0.93674 

7 

1.08633 

0.05069 

1.08633 

1.08633 

8 

1.04667 

-0.08782 

1.04667 

1.04667 

9 

0.90362 

-0.03252 

0.90362 

0.90362 

10 

0.98951 

0.10124 

0.98951 

0.98951 

Figure  13a.  Adapted  Weights  From  Amplitude  and  Phase  Nulling  in 
the  Sum  Channel 


SOLID  LINE  QUIESCENT  PATTERN 
DASHED  LINE  CANCELLATION  BEAN 


Figures  14a- 14d  show  the  results  of  placing  a  null  in  the  difference  pattern 
at  33°.  The  new  amplitude  and  phase  distributions  are  different  from  those  in 
Figure  4.  They  produce  a  null  in  the  difference  pattern  at  33°,  but  not  in  the  sum 
pattern. 


VARIABLE 

VARIABLE  WEIGHTS  TIMES 

COMPLEX  WEIGHTS 

AMPLITUDE  TAPERS 

ELEMENT 

AMPLITUDE 

PHASE 

SUM  CHAN 

DIF  CHANNEL 

1 

0.78538 

0.04267 

0.78538 

-0.78538 

2 

1.08353 

0.19367 

1.08353 

-1.08353 

3 

1.20116 

-0.07650 

1.20116 

-1.20116 

4 

0.90025 

-0.20456 

0.90025 

-0.90025 

5 

0.84762 

0.16942 

0.84762 

-0.84762 

6 

0.84762 

-0.16942 

0.84762 

0.84762 

7 

0.90025 

0.20456 

0.90025 

0.90025 

8 

1.20116 

0.07650 

1.20116 

1.20116 

9 

1.08353 

-0.19367 

1.08353 

1.08353 

10 

0.78538 

-0.04267 

0.78538 

0.78538 

Figure  14a.  Adapted  Weights  From  Amplitude  and  Phase  Nulling  in  the 
Difference  Channel 


Figure  14b.  Far  Field  Difference  Pattern  Due  to  Amplitude 
and  Phase  Nulling  in  the  Difference  Channel 


SOLID  LINE  0U1E8CENT  PATTERN 
DA8HED  LINE  CANCELLATION  BEAM 


Figure  14c.  Cancellation  Beam  From  Amplitude  and  Phase 
Nulling  in  the  Difference  Channel 


Figure  14d.  Far  Field  Sum  Pattern  Due  to  Amplitude  and 
Phase  Nulling  in  the  Difference  Channel 


These  computer  outputs  demonstrate  the  need  for  simultaneous  nulling  in  the 
sum  and  difference  channels.  Placing  a  null  in  one  channel  does  not  guarantee  that 
the  same  null  will  appear  in  the  other  channel.  When  the  simultaneous  nulling 
algorithm  was  tried,  the  results  in  Figures  15a-15e  were  obtained.  Both  antenna 
patterns  have  nulls  at  33°. 


VARIABLE  VARIABLE  WEIGHTS  TIMES 

COMPLEX  WEIGHTS  AMPLITUDE  TAPERS 


ELEMENT 

AMPLITUDE 

PHASE 

SUM  CHAN 

DIF  CHANNEL 

1 

0.77198 

-0.08627 

0.77198 

-0.77198 

2 

0.99527 

0.24139 

0.99527 

-0.99527 

3 

1.24028 

0.00000 

1.24028 

-1.24028 

4 

0.99527 

-0.24139 

0.99527 

-0.99527 

5 

0.77198 

0.08627 

0.77198 

-0.77198 

6 

0.77198 

-0.08627 

0.77198 

0.77198 

7 

0.99527 

0.24139 

0.99527 

0.99527 

8 

1.24028 

0.00000 

1.24028 

1.24028 

9 

0.99527 

-0.24139 

0.99527 

0.99527 

10 

0.77198 

0.08627 

0.77188 

0.77198 

Figure  15a.  Adapted  Weights  From  Amplitude  and  Phase  Nulling  in  the 
Sum  and  Difference  Channels  Simultaneously 


Figure  15b.  Far  Field  Sum  Pattern  Due  to  Amplitude  and 
Phase  Nulling  in  the  Sum  and  Difference  Channels 
Simultaneously 


9 


Figure  15c.  Sum  Channel  Cancellation  Beam  From 
Amplitude  and  Phase  Nulling  in  the  Sum  and  Difference 
Channels  Simultaneously 


AZIMUTH  ANGLE  IN  DEGREES 


Figure  15d.  Far  Field  Difference  Pattern  Due  to  Amplitude 
and  Phase  Nulling  in  the  Sum  and  Difference  Channels 
Simultaneously 


SOLID  LINE  QUIESCENT  PATTERN 
DASHED  LINE  CANCELLATION  BEAM 


Figure  15e.  Difference  Channel  Cancellation  Beam  From 
Amplitude  and  Phase  Nulling  in  the  Sum  and  Difference 
Channels  Simultaneously 


The  case  of  one  jammer  at  33°  was  repeated  for  phase-only  nulling.  Figures 
16a-  16d  are  the  results  of  nulling  in  the  sum  channel.  No  null  appears  in  the 
difference  pattern  at  33°.  Figures  17a-17d  show  the  results  of  difference  channel 
nulling.  Again,  the  sum  pattern  does  not  have  a  null  at  33°.  When  simultaneous 
nulling  is  used,  both  patterns  have  a  null  at  33°  (Figures  18a- 18e). 

Note  that  the  cancellation  beams  in  the  phase-only  nulling  cases  have  two  main 
lobes:  one  at  33°  and  a  replica  at  -33°.  When  added  to  the  quiescent  pattern,  a 
null  is  produced  at  33°,  but  the  sidelobe  level  at  -33°  goes  up.  The  cancellation 
lobe  at  33°  is  out  of  phase  with  the  quiescent  pattern  at  33°,  while  the  lobe  at  -33° 
is  in  phase  with  the  quiescent  pattern  at  -33°.  These  separate  beams  may  be 
mathematicallv  derived  by  substituting  Eq.  (112)  into  Eq.  (101).  The  cancellation 
beams  are  given  by 


C  BEAMS 


N 

E 


•|dn(u 


V 


M 

L 


a  sinfd  (u 
n  1  n  m 


-u  )j 
s  1 


(133) 


N  .  ,  ,  v  M 

id  (u  -  u  ) 

’  n  m  s 

a  e  7.  c  a 

,  n  m  n 

n  -  1  m  =  1 


i  E 


j  d  (u  -  u  )  j  d  (u  -  u  ) 
,  n  m  s  „  n  m  s 


/  2j 


(134) 


This  form  of  the  equation  reveals  the  two  cancellation  beams.  They  have  the  same 
amplitude,  but  conjugate  angles  and  phases. 
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ELEMENT 


PHASE  SHIFT 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


-0.18191 

0.05345 

0.16697 

-0.10013 

-0.13898 

0.13898 

0.10013 

-0.16697 

-0.05345 

0.18191 


Figure  16a.  Adapted  Weights  From 
Phase -Only  Nulling  in  the  Sum 
Channel 


Figure  16b.  Far  Field  Sum  Pattern  Due  to  PI 
Nulling  in  the  Sum  Channel 
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SOLID  LINE  QUIESCENT  PATTERN 
DASHED  LINE  CANCELLATION  BEAM 


Figure  16c.  Cancellation  Beam  From  Phase -Only  Nulling 
in  the  Sum  Channel 


Figure  16d.  Far  Field  Difference  Pattern  Due  to 
Phase -Only  Nulling  in  the  Sum  Channel 
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Figure  17a.  Adapted  Weights  From 
Phase -Only  Nulling  in  the  Sum 
Channel 


ELEMENT 


PHASE  SHIFT 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


-0.15730 

0.56268 

0.00000 

-0.56268 

0.15730 

-0.15730 

0.56268 

0.00000 

-0.56268 

0.15730 


Figure  18a.  New  Weights  Obtained  From 
Simultaneous  Sum  and  Difference  Channel 
Phase -Only  Nulling 


AZIMUTH  AN8LE  IN  DEGREES 

Figure  18b.  Far  Field  Sum  Pattern  Due  to  Simultaneous 
Sum  and  Difference  Channel  Phase -Only  Nulling 
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SOLID  LINE  QUIESCENT  PATTERN 
DASHED  LINE  CANCELLATION  BEAM 


-8S  •  “  M 

AZIMUTH  ANGLE  IN  DEGREES 


f  igure  18d.  Far  Field  Difference  Pattern  Due  to 
Simultaneous  Sum  and  Difference  Channel  Phase-Only 
Nulling 
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SOLID  LINE  QUIESCENT  PATTERN 
DA8HED  LINE  CANCELLATION  BEAM 


Figure  18e.  Difference  Channel  Cancellation  Beam  From 
Simultaneous  Sum  and  Difference  Channel  Phase -Only- 
Nulling 

The  next  set  of  computer  runs  tested  low  sidelobe  distributions  and  multiple 
jammers.  The  sum  channel  had  a  35  dB,  n  =  6  Taylor  distribution  and  the  differ¬ 
ence  channel  had  a  35  dB,  n  =  6  Bayliss  distribution.  One  jammer  was  at  23°  and 
the  other  at  50°.  First,  amplitude  and  phase  nulling  was  tried  in  the  sum  channel, 
difference  channel,  and  simultaneously.  As  in  the  previous  cases,  only  the  simul¬ 
taneous  nulling  algorithm  places  nulls  in  the  sum  and  difference  patterns  at  the 
same  time.  These  conditions  were  also  used  for  phase -only  nulling.  Figures  19 
to  25  demonstrate  the  multiple  jammer  scenario  for  a  low  sidelobe  antenna. 

A  special  case  of  interest  is  when  the  sum  and  difference  patterns  do  not  have 
the  same  sidelobe  structure.  Figures  2  6a  and  26b  show  the  results  of  simultaneous 
nulling  in  a  uniform  difference  pattern  and  a  30  dB  Chebychev  sum  pattern.  In  spite 
of  the  difference  in  sidelobe  levels  between  the  two  patterns,  each  had  a  deep  null 
at  33°.  The  resultant  Chebychev  pattern  no  longer  has  a  -30  dB  sidelobe  level,  but 
has  sidelobes  on  the  order  of  -18  dB.  From  this  example,  we  can  conclude  that 
simultaneous  nulling  works  for  sum  and  difference  patterns  with  completely  differ¬ 
ent  sidelobe  structures.  On  the  other  hand,  the  adapted  patterns  have  distorted 
mainbeams  and  sidelobe  levels. 

The  simulation  results  shown  in  this  section  demonstrate  that  the  nulling 
algorithms  work.  However,  these  results  are  only  theoretical  and  have  severe 
limitations  for  any  practical  implementations.  The  problems  and  limitations 
associated  with  simultaneous  nulling  are  discussed  in  the  next  section. 


VARIA8LE 

COMPLEX  WEIGHTS 

AMPLITUDE  TAPERS 

ELEMENT 

AMPLITUDE 

PHASE 

SUM  CHAN 

DIF  CHANNEL 

1 

1.00000 

0.00000 

0.17026 

-0.17909 

2 

1.00000 

0.00000 

0.22088 

-0.31064 

3 

1.00000 

0.00000 

0.31299 

-0.51854 

4 

1.00000 

0.00000 

0.43081 

-0.73232 

5 

1.00000 

0.00000 

0.55793 

-0.90304 

6 

1.00000 

0.00000 

0.68277 

-1.00000 

7 

1.00000 

0.00000 

0.79770 

-0.98841 

8 

1.00000 

0.00000 

0.89452 

-0.83833 

9 

1.00000 

0.00000 

0.96409 

-0.55771 

10 

1.00000 

0.00000 

1.00000 

-0.19455 

11 

1.00000 

0.00000 

1.00000 

0.19455 

12 

1.00000 

0.00000 

0.96409 

0.55771 

13 

1.00000 

0.00000 

0.89452 

0.83833 

14 

1.00000 

0.00000 

0.79770 

0.98841 

15 

1.00000 

0.00000 

0.68277 

1.00000 

16 

1.00000 

0.00000 

0.55793 

0.90304 

17 

1.00000 

0.00000 

0.43081 

0.73232 

18 

1.00000 

0.00000 

0.31299 

0.51854 

19 

1.00000 

0.00000 

0.22088 

0.31064 

20 

1.00000 

0.00000 

0.17026 

0.17909 

Figure  19c.  Quiescent  Taylor  Sum  and  Bayliss  Difference  Distributions 
(35  dB  sidelobes,  n  =  6) 


6.  LIMITATIONS  ON  SIMULTANEOUS  NULLING 

The  practical  implementation  of  the  simultaneous  nulling  algorithm  has  three 
basic  problems.  First,  the  jammer  must  be  accurately  located.  This  can  be  done 
by  using  either  a  separate  antenna  or  by  having  the  adaptive  antenna  take  time  out 
to  search  for  jammers.  The  first  idea  requires  an  extra  antenna  and  associated 
equipment.  Considering  the  cost  of  the  extra  equipment,  this  option  is  undesirable. 
The  second  idea  significantly  cuts  down  the  amount  of  time  the  radar  has  for  tracking 
and  detecting  targets.  Under  many  situations,  this  alternative  is  not  acceptable. 

In  any  case,  in  a  multiple  jammer  environment,  the  antenna  may  not  be  able  to 
accurately  locate  each  jamming  source.  Thus,  adaptive  algorithms  that  need  the 
jammer  location  in  order  to  place  the  null  are  probably  not  practical  to  implement. 

A  second  problem  with  the  algorithm  is  caused  by  hardware  tolerances  and  the 
environment.  Equipment  such  as  phase  shifters,  power  dividers,  transmission  lines, 
and  antenna  elements  create  phase  and  amplitude  errors.  Also,  environmental 
effects  like  mutual  coupling,  multipath,  and  scattering  add  unwanted  errors  to  the 
received  signals.  These  errors  appear  in  the  far  field  antenna  pattern  as  higher 
sidelobes  and  filled-in  nulls. 
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VARIABLE 
COMPLEX  WEIGHTS 


VARIABLE  WEIGHTS  TIMES 
AMPLITUDE  TAPERS 


EMENT 

AMPLITUDE 

PHASE 

SUM  CHAN 

DIF  CHANNEL 

1 

0.99812 

-0.00228 

0.16994 

-0.17875 

2 

1.00152 

0.00161 

0.22122 

-0.31111 

3 

0.99654 

-0.00481 

0.31191 

-0.51675 

4 

0.99799 

0.01049 

0.42994 

-0.73085 

5 

1.00976 

-0.00565 

0.56337 

-0.91186 

6 

0.99300 

0.00166 

0.67799 

-0.99300 

7 

1.01147 

-0.00535 

0.80685 

-0.99974 

8 

0.97940 

-0.00639 

0.87610 

-0.82106 

9 

1.00872 

0.01982 

0.97249 

-0.56257 

10 

1.00229 

-0.01262 

1.00229 

-0.19500 

11 

1.00229 

0.01262 

1.00229 

0.19500 

12 

1.00872 

-0.01982 

0.97249 

0.56257 

13 

0.97940 

0.00639 

0.87610 

0.82106 

14 

1.01147 

0.00535 

0.80685 

0.99974 

15 

0.99300 

-0.00166 

0.67799 

0.99300 

16 

1.00976 

0.00565 

0.56337 

0.91186 

17 

0.99799 

-0.01049 

0.42994 

0.73085 

18 

0.99654 

0.00481 

0.31191 

0.51675 

19 

1.00152 

-0.00161 

0.22122 

0.31111 

20 

0.99812 

0.00228 

0.16994 

0.17875 

Figure  20a.  Adapted  Weights  From  Amplitude  and  Phase  Nulling  in  the 
Sum  Channel 


Figure  20b.  Far  Field  Sum  Pattern  Due  to  Amplitude  and 
Phase  Nulling  in  the  Sum  Channel 


VARIABLE 

VARIABLE  WEIGHTS  TIMES 

COMPLEX  WEIGHTS 

AMPLITUDE  TAPERS 

ELEMENT 

AMPLITUDE 

PHASE 

SUM  CHAN 

DIF  CHANNEL 

1 

0.99858 

-0.00037 

0.17002 

-0.17883 

2 

0.99864 

0.00072 

0.22058 

-0.31022 

3 

0.99926 

0.00455 

0.31276 

-0.51816 

4 

1.00835 

0.00128 

0.43441 

-0.73843 

5 

1.00104 

-0.00849 

0.55851 

-0.90399 

6 

0.99578 

-0.00299 

0.67989 

-0.89578 

7 

0.99256 

0.00000 

0.79177 

-0.98106 

8 

0.99934 

0.00939 

0.89393 

-0.83777 

9 

1.00587 

0.00007 

0.96975 

-0.56098 

10 

1.00029 

-0.00116 

1.00029 

-0.19461 

11 

1.00029 

0.00116 

1.00029 

0.19461 

12 

1.00587 

-0.00007 

0.96975 

0.56098 

13 

0.99934 

-0.00939 

0.89393 

0.83777 

14 

0.99256 

0.00000 

0.79177 

0.98106 

15 

0.99578 

0.00299 

0.67989 

0.99578 

16 

1.00104 

0.00849 

0.55851 

0.90399 

17 

1.00835 

-0.00128 

0.43441 

0.73843 

18 

0.99926 

-0.00455 

0.31276 

0.51816 

19 

0.99864 

-0.00072 

0.22058 

0.31022 

20 

0.99858 

0.00037 

0.17002 

0.17883 

Figure  21a.  Adapted  Weights  From  Amplitude  and  Phase  Nulling  in  the 
Difference  Channel 


SOLID  LINE  QUIESCENT  PATTERN 
DASHED  LINE  CANCELLATION  BEAM 


VARIABLE 
COMPLEX  WEIGHTS 


VARIABLE  WEIGHTS  TIMES 
AMPLITUDE  TAPERS 


ELEMENT 

AMPLITUDE 

PHASE 

SUM  CHAN 

DIF  CHANNEL 

1 

0.99671 

-0.00266 

0.16970 

-0.17850 

2 

1.00016 

0.00233 

0.22092 

-0.31069 

3 

0.99577 

-0.00025 

0.31167 

-0.51635 

4 

1.00635 

0.01168 

0.43354 

-0.73697 

5 

1.01085 

-0.01405 

0.56398 

-0.91284 

6 

0.98878 

-0.00134 

0.67511 

-0.98878 

7 

1.00404 

-0.00539 

0.80092 

-0.99239 

8 

0.97868 

0.00319 

0.87545 

-0.82046 

9 

1.01459 

0.01977 

0.97815 

-0.56585 

10 

1>.00260 

-0.01378 

1.00260 

-0.19505 

11 

1.00260 

0.01378 

1.00260 

0.19505 

12 

1.01459 

-0.01977 

0.97815 

0.56585 

13 

0.97868 

-0.00319 

0.87545 

0.82046 

14 

1.00404 

0.00539 

0.80092 

0.99239 

15 

0.98878 

0.00134 

0.67511 

0.98878 

16 

1.01085 

0.01405 

0.56398 

0.91284 

17 

1.00635 

-0.01168 

0.43354 

0.73697 

18 

0.99577 

0.00025 

0.31167 

0.51635 

19 

1.00016 

-0.00233 

0.22091 

0.31069 

20 

0.99671 

0.00266 

0.16970 

0.17850 

Figure  22a.  Adapted  Weights  From  Simultaneous  Amplitude  and  Phase 
Nulling  in  the  Sum  and  Difference  Channels 


AZIMUTH  ANGLE  IN  DEGREES 


Figure  22c.  Sum  Channel  Cancellation  Beam  From 
Amplitude  and  Phase  Nulling  in  the  Sum  and  Difference 
Channels  Simultaneously 


ELEMENT 

PHASE  SHIFT 

1 

-0.00454 

2 

0.00324 

3 

-0.00959 

4 

0.02086 

5 

-0.01142 

6 

0.00334 

7 

-0.01074 

8 

-0.01251 

9 

0.03983 

10 

-0.02533 

11 

0.02533 

12 

-0.03983 

13 

0.01251 

14 

0.01074 

15 

-0.00334 

16 

0.01142 

17 

-0.01086 

18 

•  0.00959 

19 

-0.00324 

20 

0.00454 

Figure  23a.  Adapted  Weights  From 
Phase -Only  Nulling  in  the  Sum 
Channel 
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AZIMUTH  ANGLE  IN  DEGREES 


Figure  23b.  Far  Field  Sum  Pattern  Due  to  Phase -Only 
Nulling  in  the  Sum  Channel 
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Figure  23c.  Sum  Channel  Cancellation  Beam  From 
Phase-OnlyNNulling  in  the  Sum  Channel 


ELEMENT 

PHASE  SHIFT 

1 

-0.00094 

2 

0.00166 

3 

0.00898 

4 

0.00240 

5 

-0.01639 

6 

-0.00700 

7 

0.00115 

8 

0.01794 

9 

0.00055 

10 

-0.00238 

11 

0.00238 

12 

-0.00055 

13 

-0.01794 

14 

-0.00115 

15 

0.00700 

16 

0.01639 

17 

-0.00240 

18 

-0.00898 

19 

-0.00166 

20 

0.00094 

Figure  24a.  Adapted  Weights  From  Phase 
Only  Nulling  in  the  Difference  Channel 
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ELEMENT 


PHASE  SHIFT 


1 

-0.00530 

2 

0.00473 

3 

-0.00058 

4 

0.02356 

5 

-0.02857 

6 

-0.00248 

7 

-0.01074 

8 

0.00618 

9 

0.04015 

10 

-0.02796 

11 

0.02796 

12 

-0.04015 

13 

-0.00618 

14 

0.01074 

15 

0.00248 

16 

0.02857 

17 

-0.02356 

18 

0.00058 

19 

-0.00473 

20 

0.00530 

Figure  25a.  Adapted  Weights  From 
Phase-Only  Nulling  in  the  Sum  and 
Difference  Channels  Simultaneously 
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SOLID  LINE  QUIESCENT  PATTERN 
DASHED  LINE  CANCELLATION  BEAN 
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Figure  26b.  Difference  Pattern  Resulting  From 
Simultaneous  Nulling  in  a  30  dB  Chebychev  Sum 
Pattern  and  a  Uniform  Difference  Pattern 


Two  types  of  phased  array  errors  of  particular  importance  are  correlated  and 
uncorrelated  random  errors.  Uncorrelated  errors  vary  randomly  from  element  to 
element.  These  are  generally  caused  by  component  tolerances.  Correlated  errors 
result  from  groups  of  elements  in  which  each  element  has  a  related  error,  but  the 
errors  are  independent  from  group  to  group.  These  errors  slowly  vary  across  the 
aperture  and  are  due  to  uneven  heating,  bending,  edge  effects,  subarrays,  and 
so  on. 

An  antenna  pattern  consists  of  an  error  free  pattern  plus  a  pattern  due  to  errors. 

The  error  pattern  has  a  very  wide  mainbeam,  like  an  element  pattern  (cosine  shaped! 

Generally,  the  errors  are  small  and  do  not  affect  the  "high"  sidelobes  or  mainbeam. 

Thus,  close  to  the  mainbeam  the  deterministic  sidelobe  pattern  dominates  and  far 

from  the  mainbeam  the  error  pattern  dominates  (assuming  the  sidelobes  decay  going  ' 

away  from  the  mainbeam).  The  lower  the  sidelobes,  the  more  the  errors  dominate. 

The  simultaneous  nulling  algorithm  generates  a  cancellation  beam  for  a  perfect 

antenna.  This  cancellation  beam  will  not  usually  work  for  a  non-ideal  array.  The 

non-ideal  array  has  higher  average  sidelobe  levels  than  an  ideal  array.  Equation 

(135)  statistically  calculates  the  average  sidelobe  level  for  uncorrelated  random 
9 

errors. 


S  = 


(1  -  Pe)  +  ~S2  +?2Pe 
N  T)  T  Pe(  1  -7T2) 


(135) 


9.  Brookner,  E.  (n.  d. )  Antenna  array  fundamentals.  Taken  from  course  notes 

from  Microwave  Journal  Intensive  Course,  Practical  Phased-Array  Systems. 
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where 


S  =  average  far-out. sidelobe  power  level, 

N  =  number  of  elements, 

P  =  probability  of  survival  of  i^1  rading  element 
e  (calculated  from  MTBF  of  an  element), 

Tj  _  =  taper  efficiency  for  a  linear  array,  that  is, 
r  N 


whe  re 


amplitude  weight  at  element  i. 


=  mean  square  value  of  element  amplitude  error. 


_  x2  X2  -r-2 

SAMP  +  SRAD  +  ^Q 


Aamp  =  mean  square  amplifier  and  phase  shifter  fractional 
amplitude  errors. 


S  An  =  mean  square  fractional  deviation  of  radiation  element 
voltage  pattern  from  nominal  value, 

2 

2^  =  quantization  of  amplitude  weights. 


Q 

r2 


r8  =  T5 


=  mean  square  value  of  element  phase  error, 

2 


MECH  +  ^PH  +  ^POE  +  ’ 


a  =  standard  deviation  of  error.  A  gaussian  distribution  with  a  mean 
of  zero  is  assumed. 


=  mechanical  deviation  error,  or, 

MECH 


2^  ,,.,  =  (2jt/A)2  o 

MECH  MECH 


where 


A  =  wavelength  , 

2  2  ..  _  .2 
a., ......  =  (A  x  )  =  (A  v  )  , 

MEC  H  n  -  n 


Ax  =  deviation  of  element  n  in  x-direction  , 
n 

A  y  =  deviation  of  element  n  in  y-direction, 

2  2 

S’OL,  =  a  =  deviation  in  the  plane  of  the  aperture  in  the  same  units  as  A  , 

ir  ti 

2  2 

.  So  =  variance  of  a,  the  rotation  angle  error  of  the  dipole.  This 
error  is  usually  ignored  because  it  is  very  small. 


72  =  — - — irn  P  =  number  of  bits  in  phase  shifter. 
y  3X2 


In  general,  the  lower  the  sidelobes,  the  greater  the  array  amplitude  taper,  and  the 
smaller  tj  Also,  the  more  elements,  the  less  effects  the  errors  have  on  side- 
lobe  level. 

As  an  example,  consider  a  100-element  array  with  the  following  tolerances: 


<rA  =  0.043 
T)  T  =  85  percent 


CTp  =  3.58° 


~52  + 152 
(1  -C2)r]TN 


1 


0. 00185  +  0. 0039 
(1  -  0.  0039)  (0.  85)  <  100) 


0. 0000679 


=  41.  68  dB  below  peak  of  mainbeam. 


Some  other  examples  of  the  effects  of  uncorrelated  random  errors  are  found 
in  Appendix  B.  The  effects  of  uncorrelated  random  errors  on  the  simultaneous 
nulling  process  are  demonstrated  in  Figures  27  through  30.  Since  the  adaptive 
weight  changes  are  very  small  and  must  have  a  high  degree  of  accuracy,  errors 
significantly  impact  the  array's  nulling  capability.  Usually,  the  greatest  errors 
are  found  in  the  phase  shifters.  Notice  how  the  null  fills  in  as  the  number  of  bits 
decrease  in  Figures  27  and  28.  All  the  phase  shifts  are  smaller  than  the  least 
significant  bit  on  a  4-bit  phase  shifter.  The  next  figure  shows  the  pattern  after  an 
element  failure.  One  element  failure  in  a  small  low  sidelobe  array  has  devastating 
results.  Figures  30a-30c  are  cases  of  random  errors  applied  to  the  array.  These 
figures  demonstrate  the  tight  tolerances  required  to  accurately  create  a  null  in 
antenna  patterns. 
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DEGREES 


Figure  27b.  Phase-Only  Nulling  in  the  Sum  and  Difference 
Channels  Simultaneously  on  a  30  dB,  n  =  6  Taylor  Pattern 
With  9- Bit  Phase  Shifters 
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DEGREES 


Figure  27c.  Phase-Only  Nulling  in  the  Sum  and  Difference 
Channels  Simultaneously  on  a  30  dB,  n  =  6  Taylor  Pattern 
With  6-Bit  Phase  Shifters 


DEGREES 


figure  27d.  Phase-Only  Nulling  in  the  Sum  and  Difference 
Channels  Sim  janeously  on  a  30  dB,  n  =  6  Taylor  Pattern 
With  4 -Bit  Phase  Shifters 


DEGREES 


Figure  28a.  Phase-Only  Nulling  in  the  Sum  and  Difference 
Channels  Simultaneously  on  a  30  dB,  n  =  6  Bayliss  Pattern 
With  No  Errors 


Figure  28b.  Phase-Only  Nulling  in  the  Sum  and  Difference 
Channels  Simultaneously  on  a  30  dB,  n  =  6  Bayliss  Pattern 
With  8- Bit  Phase  Shifters 


DEGREES 


Figure  28c.  Phase-Only  Nulling  in  the  Sum  and  Difference 
Channels  Simultaneously  on  a  30  dB,  n  =  6  Bayliss  Pattern 
With  6-Bit  Phase  Shifters 


Figure  28d.  Phase-Only  Nulling  in  the  Sum  and  Difference 
Channels  Simultaneously  on  a  30  dB,  n  =  6  Bayliss  Pattern 
With  4 -Bit  Phase  Shifters 
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The  analysis  of  correlated  errors  is  similar  to  that  of  uncorrelated  errors. 

For  the  uncorrelated  errors  the  entire  antenna  can  be  thought  of  as  N  correlation 
regions.  In  the  analysis  of  correlated  errors,  the  antenna  is  divided  into  several 
regions.  The  errors  within  a  region  are  correlated  while  the  errors  from  region 
to  region  are  uncorrelated.  Again  we  can  think  of  the  actual  antenna  pattern  as 
being  made  up  of  several  patterns:  the  error-free  pattern,  the  pattern  due  to 
uncorrelated  errors,  and  the  pattern  due  to  correlated  errors.  Unlike  the  un¬ 
correlated  errors,  the  effects  of  the  correlated  errors  are  restricted  to  the  vicinity 
of  the  mainbeam.  Thus,  correlated  errors  are  predominant  close  to  the  mainbeam, 
and  uncorrelated  errors  are  predominant  in  the  far-out  sidelobe  regions.  Equation 

(136)  determines  the  effect  of  correlated  errors  on  the  average  sidelobe  level  of 
9 

an  array. 

s.lUJal^ilpe  /  1\  (13„ 

T)t  ng  Pe(  1  -5^)  \  VT  s  / 

where 

n  =  number  of  correlation  regions, 
s 

6  and  A  are  identical  for  all  radiating  elements  within  a  correlation  region, 

Ng  =  number  of  elements  within  correlation  region,  and 

N  =  n  N  . 
s  s 

No  computer  runs  were  made  to  demonstrate  the  effects  of  correlated  random 
errors  on  simultaneous  nulling. 

The  final  major  problem  with  the  simultaneous  nulling  arises  from  the  inherent 
limitations  of  the  phase-only  nulling  approach.  As  mentioned  earlier,  phase-only 
nulling  produces  two  equal  amplitude,  out-of-phase  cancellation  beams,  one  on  each 
side  of  the  mainbeam.  These  cancellation  begins  present  significant  problems  with 
jammers  that  are  symmetrical  about  the  mainbeam  and  jammers  in  the  mainbeam. 

Figures  3  la-3  Id  demonstrate  the  problem  with  symmetrical  jammers.  The 
phase-only  nulling  algorithm  does  not  work  for  perfectly  symmetrical  jammers. 

One  cancellation  beam  subtracts  from  the  left  side  of  the  mainbeam  and  adds  to  the 
right  side,  while  the  other  dual  cancellation  beam  tries  to  add  to  the  left  side  and 
subtract  from  the  right.  As  the  jammers  become  less  symmetrical,  the  nulling 
performance  improves.  As  a  rule  of  thumb,  phase-only  nulling  works  when  nulls 
are  not  placed  in  two  symmetrical  sidelobes. 
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Phase-only  nulling  gives  poor  results  in  the  mainbeam.  One  lobe  of  the  cancel¬ 
lation  beam  adds  to  one  side  of  the  mainbeam  and  the  other  cancellation  lobe  sub¬ 
tracts  from  the  other  side  of  the  mainbeam.  Figures  32a-32d  depict  the  results. 

The  nulling  results  from  this  section  appear  in  Table  3. 


-80  0 

DEGREES 

Figure  32b.  Cancellation  Beam  From  Phase-Only  Nulling  in 
the  Sum  and  Difference  Channels  Simultaneously  in  the  Mainbeam 
of  a  Sum  Pattern 
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Table  3.  Nulling  Results  for  Figures  27  to  32 


Far  Field  Pattern  Amplitude  (dB  below  mainbeam) 

Amplitude 

Distribution 

Jammer 

Angle 

Errors 

Far  Field 
Amplitude 

Taylor 

33° 

none 

102 

33° 

8-bit  phase  shifters 

49 

33° 

6 -bit  phase  shifters 

43 

33° 

4-bit  phase  shifters 

32 

Bayliss 

33° 

none 

101 

CO 

CO 

o 

8-bit  phase  shifters 

50 

33° 

6 -bit  phase  shifters 

40 

33° 

4-bit  phase  shifters 

31 

Taylor 

33° 

No.  2  element  failed 

26 

Taylor 

29l 

none 

112 

70° 

none 

97 

Taylor 

29° 

1  percent  amplitude, 

1°  phase. 

47 

70° 

0.  000 IX  position, 

8-bit  phase  shifters 

40 

Taylor 

29° 

2  percent  amplitude, 

4°  phase. 

16 

O 

O 

t- 

0.  001X  position, 

8-bit  phase  shifters 

42 

7.  CO  INCLUSIONS 

A  monopulse  phased  array  similar  to  the  one  in  Figure  11  requires  an  adaptive 
nulling  technique  that  places  nulls  in  its  sum  and  difference  far  field  patterns.  One 
method  of  performing  the  nulling  entails  using  one  set  of  adaptive  weights  for  the 
sum  channel  and  another  set  for  the  difference  channel  (Figure  2).  This  method 
needs  twice  as  many  adaptive  weights  and  more  computations  than  nulling  in  a 
single  channel.  Unfortunately,  nulling  in  one  channel  does  not  generate  the  same 
nulls  in  the  other  channel. 

Sections  3  and  4  outlined  techniques  for  placing  nulls  in  sum  and  difference 
patterns.  Of  particular  interest  was  simultaneous  nulling  in  the  sum  and  difference 
channels  using  one  set  of  adaptive  weights  shared  by  the  two  channels.  This  method 
offers  the  advantage  of  halving  the  number  of  adaptive  weights  needed  in  methods 
using  separate  sets  of  weights  for  each  channel.  The  disadvantages  of  the  techniques 
are  that  the  number  of  degrees  of  freedom  are  also  cut  in  half.  Since  phased  array 
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radars  usually  have  many  elements,  hence  can  have  many  degrees  of  freedom, 
reducing  the  adaptive  degrees  of  freedom  by  50  percent  will  not  appreciably  affect 
the  radar's  antijamming,  detecting,  or  tracking  functions.  However,  the  hardware 
savings  are  substantial. 

The  problem  of  practical  implementation  of  the  simultaneous  nulling  technique 
needs  further  investigation.  The  technique  described  does  not  employ  feedback 
to  adjust  the  nulls  for  the  nonideal  antenna  pattern.  The  technique  most  be  modified 
to  incorporate  a  feedback-control  system  to  adjust  the  adaptive  weights.  One  way 
to  change  the  algorithm  id  to  adaptively  adjust  the  cancellation  beams  until  they  are 
at  the  same  height  as  the  nonideal  pattern,  rather  than  the  ideal  pattern.  Now, 
when  the  cancellation  beams  and  nonideal  pattern  are  added,  the  resultant  pattern 
will  have  a  null  in  the  desired  direction. 

The  phase  shifts  required  to  place  nulls  in  the  antenna  pattern  are  quite  small 
and  prone  to  component  errors.  Phase  shifter  quantization  is  probably  the  major 
factor  limiting  the  accuracy  of  nulling.  Building  digital  phase  shifters  with  more 
than  6  or  8  bit  control  lines  is  not  a  reasonable  undertaking  at  this  time.  They 
would  be  too  expensive  to  be  practical  and  the  last  few  bits  would  probably  be  in  the 
noise  level. 

An  alternative  to  the  many-bit  phase  shifter  is  two  separate  sets  of  phase 
shifters.  One  set  would  make  coarse  adjustments  for  steering  the  main  beam.  They 
would  have  4  bits  or  22.  5°  accuracy.  The  other  set  would  be  phase  shifters  which 
could  be  accurately  set  between  plus  and  minus  2°.  These  phase  shifters  would 
provide  the  small,  accurate  phase  adjustments  needed  to  place  the  nulls  in  the 
antenna  patterns. 

Besides  the  topic  of  practical  implementation,  other  areas  in  simultaneous 
nulling  need  further  investigation.  For  instance,  partial  adaptive  nulling  seems  to 
be  an  appealing  way  to  reduce  the  number  of  adaptive  weights,  while  maintaining  the 
nulling  capability.  This  technique  will  probably  work  well  theoretically.  Another 
concern  is  mainbeam  nulling.  What  are  the  limitations  and  how  is  the  monopulse 
ratio  affected?  Finally,  can  simultaneous  nulling  be  used  without  cutting  the  degrees 
of  freedom  in  half?  Perhaps  there  are  sum  and  difference  pattern  distributions 
such  that  when  a  null  is  placed  in  one  pattern,  the  null  also  appears  in  the  other 
pattern.  These  are  just  a  few  of  the  areas  of  theoretical  research  that  should  be 
pursued. 

This  report  has  described  the  importance  of  nulling  in  both  the  sum  and  differ¬ 
ence  patterns  of  a  monopulse  antenna.  Nulling  in  either  the  sum  or  difference 
channels  will  not  place  the  same  null  in  both  far  field  patterns.  It  is  theoretically 
possible  to  simultaneously  place  nulls  in  both  far  field  patterns  by  adjusting  one  set 
of  adaptive  weights  shared  by  the  two  channels.  If  this  technique  can  be  used  on 
phased  .irray  antennas,  there  will  be  substantial  savings  in  hardware.  First,  a 
practical  way  to  implement  the  technique  on  actual  phased  array  antennas  is  needed. 
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Appendix  A 


Computer  System 

The  low  sidelobe  synthesis  techniques  and  nulling  algorithms  were  programmed 
in  BASIC  on  the  Tektronix  4052  computer.  A  diagram  of  the  computer  set  up  is 
shown  in  Figure  Al.  The  Tektronix  4052  has  64K  of  memory  and  tape  storage 
capability.  To  the  left  of  the  computer  is  the  Tektronix  4663  plotter.  It  drew  all 
the  plots  used  in  this  report.  On  the  right  is  a  Tektronix  4631  hard  copy  unit  and 
Tektronix  4641  printer. 

Figure  A2  shows  a  block  diagram  of  the  computer  algorithm.  First,  data  such 
as  the  number  of  elements  and  element  spacing,  are  entered  on  the  keyboard.  Next, 
the  program  calculates  the  amplitude  weights  for  either  a  uniform,  Taylor,  or 
Chebychev  sum  distribution  and  for  either  a  uniform  or  Bayliss  difference  distribu¬ 
tion.  This  part  of  the  program  models  the  low  sidelobe  synthesis  techniques  des¬ 
cribed  in  Sect. on  2.  Multiplying  the  zero  factors  together  to  obtain  the  polynomial 
p(s)  is  very  inefficient.  Only  linear  array  distributions  with  up  to  30  elements  can 
be  derived  without  taking  excessive  computer  time. 

After  calculating  the  desired  array  distribution,  the  program  enters  into  the 
nulling  routines.  The  routine  offers  six  different  choices: 

(1)  Amplitude  and  phase  nulling  in  sum  channel, 

(2)  Amplitude  and  phase  nulling  in  difference  channel, 

(3)  Amplitude  and  phase  nulling  simultaneously, 

(4)  Phase  only  nulling  in  sum  channel, 

(5)  Phase  only  nulling  in  difference  channel, 

(6)  Phase  only  nulling  simultaneously. 


The  routine  changes  the  quiescent  sum  and  difference  channels  for  one  of  the 
above  choices.  Both  channels  are  assumed  to  share  one  set  of  adaptive  weights. 

In  cases  1  and  4  the  adaptive  weights  are  adjusted  to  place  nulls  in  the  sum  channel. 
The  difference  channel  receives  the  exact  same  weight  change.  For  cases  2  and  5 
the  reverse  holds  true.  In  cases  3  and  6  the  adaptive  weights  are  adjusted  to  yield 
nulls  simultaneously  in  the  sum  and  difference  channels.  When  the  nulling  routine 
is  completed,  the  program  stores  the  quiescent,  adapted,  and  cancellation  beam 
weights  on  tape.  Next,  the  program  calculates  the  far  field  pattern  of  any  of  the 
array  weights  stored  on  tape.  The  user  interactively  specifies  the  number  of  plots 
on  the  graph,  the  number  of  far  field  transform  points  on  each  plot,  the  angular 
extent  of  the  graph  and  the  normalization  value.  In  this  way,  the  output  can  be 
modified  to  obtain  the  necessary  details.  Finally,  the  computer  draws  the  graph 
on  the  CRT  and/or  the  plotter. 

I  wrote  all  parts  of  the  program  except  for  the  complex  inversion  routine.  That 
routine  came  from  a  Tektronix  library  tape. 
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Figure  A  1.  Block  Diagram  of  the  Computer  Layout 
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Figure?  A2.  Block  Diagram  of  the  Computer  Program 


Appendix  B 


Effects  of  Errors  on  Sdotobo  Levels 

The  following  pages  contain  diagrams  of  RMS  sidelobe  level  vs  number  of  an¬ 
tenna  elements  for  errors  due  to 

(1)  Random  phase  variations, 

(2)  Random  amplitude  variations, 

(3)  Phase  shifter  quantization, 

(4)  Random  variations  of  element  position, 

(5)  Element  Failures. 

These  uncorrelated  random  errors  are  the  inherent  limitations  to  adaptive 
nulling.  Better  nulling  requires  better  component  tolerances,  especially  in  low 
sidelobe  antennas. 

The  errors  in  diagrams  B1  to  B15  were  calculated  from  Eq.  (135).  A  normal 
distribution  was  assumed  for  the  phase,  amplitude,  and  position  errors.  These 
errors  have  standard  deviations  of  O'  a.,  and  <7  respectively.  Element  failures 

x'rl  A  p  ^ 

are  represented  by  Pe«  the  probability  of  survival  of  the  i  radiating  element.  The 
effects  of  errors  (including  element  failures)  go  down  as  the  number  of  elements, 
gain,  and  aperture  efficiency  go  up. 


97 


NUMBER  OF  ANTENNA  ELEMENTS 


Figure  Bl.  RMS  Sidelobes  Due  to  Random  Phase 
Errors  (jj  =  1.  0) 
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Figure  B5.  RMS  Sidelobes  Due  to  Random  Amplitude 
Errors  (tjt  =  0.  8) 
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Figure  B6.  RMS  Sidelobes  Due  to  Random  Amplitude 
Errors  (tj-  *  0.  6) 
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Figure  B9.  RMS  Sideiobes  Due  to  Quantization 
Errors  *  0.  6) 
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Figure  Bll.  RMS  Sidelobes  Due  to  Element  Position 
Errors  (tj^  =  0.  8) 


Figure  B12.  RMS  Sidelobes  Due  to  Element  Position 
Errors  (fyp  =  0.  6) 
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